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ABSTRACT 

We use a sample of 47 homogeneous and high sensitivity CO images taken from the Nobeyama and 
BIMA surveys to demonstrate that, contrary to common belief, a significant number (~ 40%) of HI- 
deficient nearby spiral galaxies are also depleted in molecular hydrogen. While Hl-deficiency by itself 
is not a sufficient condition for molecular gas depletion, we find that H2 reduction is associated with 
the removal of HI inside the galaxy optical disk. Those Hl-deficient galaxies with normal H2 content 
have lost HI mainly from outside their optical disks, where the H2 content is low in all galaxies. This 
finding is consistent with theoretical models in which the molecular fraction in a galaxy is determined 
primarily by its gas column density. Our result is supported by indirect evidence that molecular 
deficient galaxies form stars at a lower rate or have dimmer far infrared fluxes than gas rich galaxies, 
as expected if the star formation rate is determined by the molecular hydrogen content. Our result is 
consistent with a scenario in which, when the atomic gas column density is lowered inside the optical 
disk below the critical value required to form molecular hydrogen and stars, spirals become quiescent 
and passive evolving systems. We speculate that this process would act on the time-scale set by the 
gas depletion rate and might be a first step for the transition between the blue and red sequence 
observed in the color-magnitude diagram. 

Subject headings: Galaxies: evolution - galaxies: ISM - galaxies: spiral - ISM: molecules - stars: 
formation 



1. INTRODUCTION 

A complete and coherent description of the detailed 
connection between a galaxy's gaseous phase and its star 
formation is still lacking. This poses a significant lim- 
itation for studying galaxy evolution, even in the local 
universe. From an observational point of view, two dif- 
ferent approaches have been adopted to address this fun- 

' damental topic. The first takes advantage of high spatial 
resolution images to compare directly the molecular and 
atomic phases of the interstellar medium (ISM) with re- 
cent star formation activities. The results (W ong & Blita 

; 12002*: 'Ken nicutt et al.l 120071 : iFumagalh fc Cavazzi.. 20081 

. [Bigiel et al.l l2008[ ) suggest that stars form preferentially 
in molecular regions, as traced by diffuse CO emission. 

j The second approach is more indirect, but relies on 
estimators for the star formation rate and atomic gas 
that can be easily collected for large samples. These 
projects statistically study the effects of the atomic 
gas removal on star formation in cluster galaxies com- 
pared with field galaxies. Results (Kennicutt & Kent 
1983; Gavazzi et alj 120021 : iKoopmann fc KennevI 120041 : 
Gallazzi et al.. ,20091 Gavazzi et. al, in prep.) show a 
significant correlation between the atomic gas and the 
star formation: the latter is quenched in those galax- 
ies which are HI deficient due to environmental effects. 
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Although this result ca n be interpreted in terms of HI 
feeding star formation (jBoselli et al.ll200lf) . it is possi- 
ble that molecular hydrogen (H2) depletion also plays a 
critical role. 

However, observations to date suggest that molec- 
ular hydrogen, which dominates the mass of the in- 
terstellar medium in the center of spiral galaxies, is 
not strongly affected by environmental perturbations. 
Galaxies in rich clusters appear to have on average the 
same molecular content as fi eld g a laxies (fSt ark et saL] 
1986; Kcnnev & Young' (19861 11981 ICasoh et al.1 119p 
Boselli et al. 1997. 2002i) . This result is consistent with 
a simple theoretical argument based on the assumption 
that ram-pressure stripping is the dominant environ- 
mental process, a fact that, although not yet proved, 
is likely to be the correct interpretation for HI obser- 
vat ions in the Virgo and Coma clusters (see a review 
bv lBoselli fc Gavazzill2006[ l. According to this picture, 
a galaxy which travels at a velocity Uroi relative to 
the dense intergalactic medium (IGM) or intracluster 
medium (ICM) characterised by a density pioM looses 
its ISM at the escape velocity Vesc and density pisM only 
if 

PlGMt^rcl > PlSM^Lc ■ (1) 

This condition is not easily satisfied for gas in molecular 
hydrogen phase, which has densities up to 10^ cm""^, 
and is usually strongly bound to the galaxy due to its 
proximity to the galactic center. 

This theoretical argument seems to agree with observa- 
tions of the CO abundance in members of the Virgo and 
Coma cl usters. Using a sa mple composed of 47 Virgo 
galaxies, ^Star k et al.l ( 19861 ) showed that the molecular 
fraction CO/HI increases towards the center of the clus- 
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TABLE 1 
Properties of our galaxy sample. 
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R.A. 


Dec. 


Dist. 
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-0.24 
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59.60 


3.4 
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10.1 
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-0.17 
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iNCjrO Uzoo 


11.89 
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3.0 
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9.7 
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35 


13 * 

r> 
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14.4 
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50 
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143.04 
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6.3 
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10.4 
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0.33 


18 
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154.57 


41.42 


8.7 


Sc 


10.2 


222 


0.11 


0.32 


- 


N* 


■\T r~\ r~\ o o er 1 


160.99 


11.70 


8.1 


Sb 


10.4 
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0.34 


0.60 


16 


N* 


TVT /"I o cr A /I 


165.80 


27.97 


26.5 


Sab 


10.7 


73 


0.10 


0.37 


35 
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AT O C O 1 


166.45 


-0.03 


7.2 


Sb 


10.6 


249 


0.19 


0.06 


20 


N* 


JNCjrO OOZY 


170.06 


12.99 


6.6 


Sb 


10.5 


306 


-0.11 


0.81 


19 


N* 


1n<jL/ oboi 


170.26 


53.17 


21.6 


be 


10.7 


111 


0.01 


-0.03 


34 


IN 




178.20 


44.12 


17.0 


Sc 


10.5 
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-0.06 


-0.24 


26 


B 


AT /I n c; 1 
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44.53 


17.0 


Sb 


10.6 


146 


-0.03 


0.26 


32 


N 
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181.60 


52.71 


17.0 


Sb 


10.5 


88 


-0.02 


0.47 


26 


N 


AT /lino 
JNCjrO 4192 


183.45 


14.90 


17.0 


Sb 


11.0 


293 


0.59 


0.36 


9 


N 


AT /I o 1 o 


183.91 


13.90 


16.8 


Sc 


10.3 


108 


0.17 


0.57 


20 


N 


AT /I O CT /I 

JNCjrO 4254 


184.71 


14.42 


17.0 


Sc 


10.9 


184 


-0.20 


0.18 


32 


N 


iNLxO 42oo 


184.74 


47.30 


6.8 


bb 


10.7 


545 


0.44 


0.25 


15 


15 
r> 


JNCjO 4jUo 


185.48 


4.47 


17.0 


Sc 


11.0 


197 


-0.08 


0.13 


36 


N 


AT /I o o 1 
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15.82 
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Sc 


11.1 


273 


0.05 


0.53 


18 


N* 


AT/^ /-I yi fir> 
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186.53 


13.11 


17.0 


be 


10.4 


118 


0.02 


0.73 


12 


IN 


AT/^ A A 


186.61 


31.22 


17.0 


be 


10.9 


108 


-0.09 


-0.15 


23 


IN 


AT /I /I 1 n 

JNCjrO 4419 


186.74 


15.05 


17.0 


Sa 


10.6 


105 


0.40 


1.04 


7 


N* 


iNCjrO 4oUl 


188.00 


14.42 


17.0 


bbc 


11.2 


216 


0. 14 


0.51 


6 
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IN 


■\T/^/-i /I croc: 


188.58 


8.20 


17.0 


be 


10.8 


249 


0.02 


0.35 


17 


IN 




188.61 


2.19 


17.0 


Sc 


10.7 


216 


0.18 


0.29 


20 


N* 


AT /I cr /I o 


188.86 


14.50 


17.0 


Sb 


10.9 


180 


0.84 


0.70 


17 


N 


INCjrO 45d9 


189.21 


13.16 


17.0 


bab 


11.0 


321 


0.23 


1.13 


2 


IN 


JNCjrO 4579 


189.43 


11.82 


17.0 


Sab 


11.1 


188 


0.62 


0.71 


4 


N 


AT Ant^ A 
JNCjrO 4o54 


190.99 


13.13 


17.0 


Sc 


10.7 


149 


-0.02 


0.17 


31 


N* 


INijrO 4ooy 


191.94 


13.76 


17.0 


be 


10.5 


175 


0. 14 


1.06 


14 


IN 


IN 4 / OD 


192.72 


41.12 


4.3 


bab 


10.4 


232 


0.78 


0.68 


10 


IN 


JNCjC 4o2o 


194.18 


21.68 


4.1 


Sab 


10.3 


314 


0.70 


0.86 


9 


B* 


JNOrO 5(JU5 


197.73 


37.06 


21.3 


Sb 


11.3 


143 


0.28 


0.47 


7 


B 


iNCjrO OUoo 


198.36 


36.59 


18.7 


be 


11.0 


293 


-0.16 


0.19 


20 


15 
r> 


NGC 5055 


198.96 


42.03 


7.2 


Sbc 


10.7 


352 


0.04 


0.12 


20 


N* 


NGC 5194 


202.47 


47.19 


7.7 


Sbc 


10.8 


232 


-0.36 


0.12 


28 


N* 


NGC 5236 


204.25 


-29.86 


4.7 


Sc 


10.7 


464 


0.23 


-0.49 


33 


N* 


NGC 5247 


204.51 


-17.88 


22.2 


Sb 


10.9 


161 


0.06 


0.05 




N 


NGC 5248 


204.38 


8.88 


22.7 


Sb 


11.1 


122 


-0.10 


-0.25 


26 


N 


NGC 5457 


210.80 


54.35 


5.4 


Sc 


10.5 


719 


0.06 


0.28 


31 


N* 


NGC 6217 


248.16 


78.20 


23.9 


Sbc 


10.5 


67 


0.13 


-0.45 


34 


N 


NGC 6946 


308.72 


60.15 


5.5 


Sc 


10.4 


344 


-0.28 


-0.45 


34 


N* 


NGC 6951 


309.31 


66.10 


20.3 


Sb 


11.1 


97 


0.24 


0.04 


16 


N 


NGC 7331 


339.27 


34.42 


14.3 


Sbc 


11.1 


273 


0.14 


-0.04 




B* 


UGC 2855 


57.09 


70.13 


20.3 


Sc 


10.9 


132 


-0.34 


-0.13 


34 


N 



References. — References are discussed through 
^ N: Nobeyama survey (Kuno et al."200,7'); B: BIMA 
high— resolution study are identified by an asterisk. 

ter, identified by the position of M87. This is consis- 
tent with the removal of atomic hydrogen only from in- 
ner cluster galaxies, without any corresponding change 
in the molecular hydrogen content. With independent 
CO measurement s, the same c onclusion was reached by 
IKennev fc YouneJ (jl986l Il989f ). who studied molecular 
hydrogen deficiency in datasets of 23 and 41 Virgo galax- 
ies. Their analysis shows that the molecular content 
has not responded to the HI removal on time scales of 
10^ yr, i.e. it is unpe rturbed with i n the Virgo cross- 
ing time. Furthermore, ICasoh et al.l (|199in pointed out 
that even in the Coma cluster, where the ram-pressure 
is stronger than in Virgo, there is no correlation be- 
tween the H2 mass and the atomic hydrogen content 
of individ ual spirals. Surpr is ingly, using the data col- 
lected by IKennev fc Youn^ (|1989f ) with updated dis- 



out the text. 

survey IIHelfer et al.ll2003) . Galaxies included in the 

tances, IRengaraian fc Ivengail (|1992f ) found from a com- 
parison of normalized H2 masses that some deficiency of 
molecular hydrogen holds among Virgo galaxies, as com- 
pared to the field. This re s ult is in con trast with a series 
of papers bv lBoselh etHI (|T997ll200^ . who did not find 
any correlation between the HI and the H2-deficiency 
on a large sample of ^ 250 Virgo and Coma galaxies. 
Except some extreme cases of galaxies w hich have their 
mole cular gas depleted and displaced (e.g. lVollmer et all 
l2008f) . there is no strong evidence that cluster member 
galaxies have statistically less molecular gas than field 
galaxies. 

This seems in contradiction with a coherent picture 
of the star formation in spirals. If stars really do form 
primarily in molecular gas, one might expect the quench- 
ing in the star formation activity to be connected with 
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a molecular gas deficiency in HI p oor galaxies. This hy- 
pothe sis is partially supported by iFumagalli fc Gavazzil 
((200l) . who, using a sma ll sample of 28 im ages from 
the Nobeyama CO survey (iKuno et al.ll2007f l. find some 
hints of a possible depletion in the molecular content 
(see their Fig. 4 and 6). The aim of this paper is test for 
H2-deficiency by using homo geneous and high sensitiv- 
ity data from the Nobeyama (jKuno et aIll2007D . BIMA 
(iHelfer et al.|[20 03l. VI VA (C hung at al., in prep.) and 
THINGS (|miter et al.i 200g) surveys. At first, we com- 
pare the integrated molecular and atomic gas mass; then, 
we study the distribution of the g clS 3tS cl function of the 
galaxy radius and finally we compare both the atomic 
and hydrogen gas content with indicators for the star for- 
mation. While previous studies selected galaxies mainly 
from the Virgo and Coma clusters, our sample includes 
also objects from nearby groups. In these diverse envi- 
ronments it is less obvious that ram-pressure is the most 
efficient process in removing the atomic gas. Therefore, 
in our analysis we will focus only on the effects produced 
by the HI removal, without addressing to the details of 
its origin that is likely connected with a generic environ- 
mental perturbations or galaxy feedback. 

The paper is organized as follows: the sample and the 
data reduction are discussed in Sec. [U the analysis is 
given in Sec. [3] and the discussion in Sec. 01 In par- 
ticular, we focus on why previous surveys did not detect 
the H2-deficiency ([23]), we explain the process responsi- 
ble for the molecular gas depletion (i i4.2p and we discuss 
the implications of H2-deficiency on the star formation 
fi j4.3p . The conclusions are in Sec. El 

2. DATA REDUCTION 
2.1. The sample 

Our sample is selected based on the criterion of high 
sensitivity CO images and homogeneity to minimize the 
scatter in CO measurements. Because previous CO sur- 
veys show a relative flux dispersion which severely lim- 
its the H2 deficiency analy sis fsee |j 4. ip. we limit to 
the Nobeyama CO survey ( Kuno et al.l 120071 ) and the 
BIMA survey (|Helfer et all I2003D which are shown to 
agree to within 0.4 dex. 47 massive spiral galaxies in 
the H-band luminosity range 10^° < Lh/Lq < 10^^ 
meet our selection criterion, 40 from the Nobeyama and 
7 from BIMA. Our sample includes only spiral galaxies 
since early-type (E /dE) galaxies are known to be ex- 
tremely gas poor ((di Serego Ahghieri et al.ll2007f) . while 
late-typ e galaxies (Sa-Sc. BC D or Irr/dirr) are typically 
gas rich (jCavazzi et al.l[2b08[ ). Table [T] summarizes some 
observational properties for our sample. Since we are 
dealing mainly with the full Nobeyama CO atlas, which 
has been selected using a far infrared (FIR) emission 
criterion (i^ioo/^m > 10 Jy), our sample may be biased 
against CO poor galaxies. 

The Nobeyama maps are obtained with a 45m single- 
dish radio telescope that provides high sensitivity ('^ 0.5 
Mq pc^) flux measurements for the entire galaxy disk. 
The spatial resolution is about 15", which corresponds to 
an average physical scale of ~ 1 kpc at the mean distance 
of our sample 15 Mpc). In contrast, despite its higher 
spatial resolution (6" — 10", corresponding to 450-700 
pc) , the interferometric BIMA survey has systematically 
poorer sensitivity to extended flux because interferome- 



ters cannot detect flux from regions larger than the size- 
scale set by the minimum distance between the antennas. 
One can recover the total flux by using single dish obser- 
vations that provide the zero spacing element. We have 
limited our galaxy sample to the subset with single-dish 
data. Even for these galaxies, the r econstruction algo- 
rithm is biased against large scales (jHelfer et al.ll2003L 
app. A); however, a comparison between Nobeyama and 
BIMA data indicates that the missed flux is negligible 
compared to the total flux, i.e. the BIMA survey provides 
a reliable estimate of the integrated molecular hydrogen 
mass. When a galaxy is imaged in both the surveys, we 
use the Nobeyama data. 

For our analysis we will also compare the spatial dis- 
tributions of molecular and atomic hydrogen; therefore 
we colle ct spatially resolve d HI images from the THINGS 
survey (jWalter et al.ll200 8^ and the VIVA survey (Chung 
et al., in prep.). Since spatially resolved images are not 
available for all the galaxies in our sample, while study- 
ing surface density proflles we restrict our analysis to a 
subsample of 20 galaxies (see Table [TJ. 

2.2. A multifrequency analysis 

The comparison between the molecular distribution, 
the atomic hydrogen content and the star formation ac- 
tivity requires some preliminary transformations from 
observed to physical quantities. In this section we de- 
scribe this issue. 

2.2.1. The molecular hydrogen 

The H2 molecule has no permanent dipole so we must 
infer its abundance indirectly by assuming that the con- 
ditions of the ISM where the H2 forms are similar to 
those in which other molecules are synthesized. A good 
tracer for H2 is CO because the intensities of its rota- 
tional transitions are correlated to the molecular column 
density. This is based on a simple argument which re- 
lates the velocity inferred from th e line profile with the 
virial velocity of mo lecular clouds (jPickman et aLlfigSGl : 
iSolomon et"allll987f l. 

The Nobeyama and BIMA surveys image the CO 
J : 1 — transition at 115 GHz, which is excited above 
5.5 K in cool gas at moderate density. We infer the 
molecular hydrogen column density by adopting a con- 
version factor X (cm~^ (K km/s)^^), which is a func- 
tion of temperature, density, UV radiation field, metal- 
licity and even of the shape of the molecular clouds ( Wal| 
l2007HBosem et al.ll2002lll997l app. A). Because it is im- 
possible to model X theoretically for actual galaxies, a 
constant conversion factor derived using observations in 
our G alaxy has often been assumed (e.g. ISolomon et al.l 
119871 ). The Galactic factor likely gives underestimates 
for the molecular hydrogen content in low mass galaxies, 
because empirical ca fibrations show tha t X decreases for 
increasing mass (e.g. lBoselli et al.ll2002[ ). Since we expect 
the CO-to-H2 conversion factor to depend primarily on 
the dust and metallicity that regulate the formation and 
photodissociation of H2 and the CO heating, we choose 
an individual X valu e for each gal a xy fol lowing the em- 
pirical calibrations of iBoselli et al.l (|2002f ). This is 

logA: = -0.381ogLH + 24.23, (2) 

where Lh is the H-band luminosity in L©. We retrieve 
Lh from GOLDMine (Gavazzi et al.,200al and SIMBAD. 
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This calibration accounts for the dependence of X on 
the metal abundance because the H- band luminosity cor - 
relates with galaxy metallicity fe.g. iBoselli et alll2002f ). 
The systematic uncertainty on the X factor derived using 
eq. ([2]) ranges from 0.65 to 0.70 dex for the luminosity in- 
terval 10^*^ — 10^^ Lq. Other unknown dependencies can 
increase this error. The final conversion for the molecular 
hydrogen surface density is: 



(3) 



where mp is the proton mass in Mq and _Fco is the 
CO flux in K km/s. This conversion does not include 
the mass of helium, but its contribution may be in- 
cluded trivially by multiplying the derived surface den- 
sity by 1.36. The integrated CO flux {Fco in Jy km/s) 
is converted into a total H2 mass (Mhj in Mq) using 
(IHelfer et al. l [2003l) : 

Mh^ = 3.92 X IQ-^'^XFcoD^ , (4) 

whe re D i s the g alaxy distance in Mpc from GOLDMine 



wiie re L) i s tnega 
and|™i3 (fT99l . 



2.2.2. The atomic hydrogen 

We estimate the atomic hydrogen content via 21 cm 
emission line. Single-dish HI s urveys, inclu ding the AL- 
FALFA (|Giovanelli et al.l[2005h or HIPASS (|Mever et all 
[200I surveys, do provide the HI integrated flux for all 
the galaxies in our sample. Collecting different HI data 
for a same galaxy, we notice a signiflcant dispersion in the 
measured HI flux, even up to a factor of 2. Because we 
are interested in a relative comparison between galaxies 
rather than an accurate determination of the absolute HI 
mass, we re trieve all the HI flu xes from the HyperLeda 
database (Paturel et al.l l2003aD . Even if more recent 
measurements have better sensitivity and provide a more 
accurate determination of the total HI ab undance, this 
approa ch gives an homogeneous dataset (jPaturel et al.l 
l2003bf) . As discussed in Appendix \M a particular choice 
for the observed flux does n_ot affect our analysis. 

The integrated HI flux (Fhi in Jy km/s) is converted 
into a total mass (Mhi in Mq) using 

Mm = 2.36 X IO^Fhi^*^ , (5) 

whe re D i s the galaxy distance in Mpc from GOLDMine 
and iTullvl (1994) . The HI surface density profile is com- 
puted from THINGS and VIVA images according to the 
conversion 



Eh/ [Mq cm 2] = mpiVHi 



(6) 



where iVm is the column density in cm ^ derived from 
the flux J ^HT (Jy/beam km/s) using the following calibra- 
tion (e.g. iWalter 61111 12008): 



A^Hi[10^^cm^] = 110.4 X 10-^ 



HI 



n 



(7) 



with B the beam FWHM in arcsec. As with H2, these 
surface densities do not include the mass of helium. 

2.2.3. The star formation activity 

In the ISM of galaxies, the Ha emission line is caused 
by the recombination of atomic hydrogen which has been 
ionized by the UV radiation of O and B stars. Therefore 



Ha is a tracer of recent, massive {M > 10 Mq) star for- 
mation activity on time scales oft < 4x 10^ yr ( Kenniciittl 
119981 ). For all but three o bjects, we retrieve from 
the literature ( GOL Dmine: | Kennicutt fc Kent! 119831 : 
iJames et aL|[2004l: iKennicutt et al.ll2008f) the Ha squiva- 
lent width (E.W.), a tr acer of the speciflc star formation 
rate (|Kennicuttlll983b[ ). 

Since the narrow band fllter used to image the Ha 
emission line is large enough to collect the flux from 
the [Nil] line, the measured E.W. overestimates the 
Ha emission. Unfortunately, the ratio [Nil]/ Ha is not 
available for all of the galaxies in our sample and we 
can apply only t he average c orrec tion [Nil]/ Ha = 0.53, 
as suggested by IKennicuttI (jl992D . However, this value 
underestimates the [Nil] emission from the bulges of 
galaxies that host active galactic nuclei (AGN), where 
the nuclear [Nil] can be even stronger than the Ha. 
In addition, even if the nuclear emission were corrected 
using a suitable factor, it would not be possible to disen- 
tangle between the Ha flux due to star formation from 
that due to AGN ionization. Our sample, composed of 
massive s pirals , is lik ely to be affected by this contami- 
nation ([Decarli et al.li2007il . 

Another effect to consider is dust absorption. In prin- 
ciple one can estimate the extinction by measuring the 
Ha/ Hp or the Ha/Br^ fluxes (case B recombination). 
These values, however, are unknown for the majority of 
our galaxies. One could apply a stand ard dust ex tinc- 
tion (e.g. A = 1.1 mag. IKennicutt fc Ken^ I1983D but 
this would not change the relative distribution of the 
measurements. Instead, we simply keep the measured 
Ha -f [Nil] E.W. as an approximation for the specific 
star formation rate. The reader should keep in mind 
that a large uncertainty is associated with these measure- 
ments, in addition to the error in the flux determination. 

2.3. Global fluxes and surface brightness profiles 

We compute total CO fluxes from the Nobeyama 
images using the IRAF task qphot by integrating 
the emission in an aperture which encloses the entire 
galaxy. Fluxes from the BIMA survey are retrieved from 
IHelfer et al.l (|2003^, who list values computed within a 
square region. We retrieve the HI flux from HyperLeda, 
as previously discussed. We compute the surface bright- 
ness profiles for the CO and HI using an IDL procedure 
that averages the flux into elliptical annuli characterised 
by three fixed parameters (the galaxy center, the eccen- 
tricity and the position angle) ; we adopt a step along the 
major axis of ~ 5 arcsec, roughly 1/3 of the spatial reso- 
lution of the CO images from the Nobeyama survey. All 
the profiles are corrected for the projection effect using 
the galaxy inclination (HyperLeda) on the plane of the 
sky. For each galaxy we finally obtain the HI and H2 
surface density profiles Yihi and E/f^ iii -M© pc~^. 

3. ANALYSIS 

To test for molecular hydrogen depletion in our 
dataset, we compare the H2 mass in a sample of iso- 
lated galaxies with that in perturbed galaxies. Our iso- 
lation criterion is defined by the atomic hydrogen con- 
tent. Because the edges of the HI disk are weakly 
bounded, any perturbation can easily remove atomic 
gas from a galaxy, lowering its HI mass. That the 
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Fig. 1. — a - Histogram of the Hl-deficiency values in our sample, separated in isolated and perturbed galaxies at de/jji = 0.3. A positive 
value of defm indicates that a galaxy has an HI mass lower than expected on average according to its size and morphology, b - Histogram 
of the H2-deficiency values for isolated galaxies only. The Gaussian fit shows that a = 0.2. 



Hl-deficie ncy should trace environmen tal effects is sup- 
ported by IGiovanelli fc HavnesI (|1985f ). who show that 
Hl-deficiency increases toward the center of a cluster. In 
these regions the IGM is dense, the cluster potential well 
is deep and the number of galaxies is high so that en- 
vironmental processes are very effective. We expect HI 
deficiency to be a useful measure of disturbance even if 
environment is not the dominant means of gas removal. 
For example, if star formation feedback is the primary 
cause of gas loss, this too will preferentially remove the 
low-density, loosely bound HI. 

Formally, the atomic gas depletion is quantified 
using defm, the deficien cy parameter, defined by 
IHavnes fc Giovanellil (jl984f ) as the logarithmic difference 
between the expected HI mass in isolated galaxies and 
the observed value: 



Only 4 galaxies lie at de/ni < —0.3, indicating that our 
reference sample is not dominated by galaxies with HI 
mass above the average. The remaining 22 galaxies have 
defm > 0.3, spanning the entire range of moderate and 
high deficiencies (0.3 — 1.2 dex). 

In a similar way, we quantify a possible reduction in 
the molecular hy drogen content using the H2-deficiency 
parameter de/n^ (IBoseUi et al.lll997l l2002l ): 



defn^ = logMn^ 



l0gMH2,obs 



(9) 



defm = log A/hi, cxp - logA/Hi,obs 



(8) 



The reference mass Mhi,cxp is computed as a function 
of the optical radius at the 25th mag/arcsec^ in B band 
(retrieved from NED and HypcrLcd a) us i ng th e coeffi- 
cients recently updated by Solancs et al.l (jl996f ). which 
are weakly dependent on the morphological type. Sev- 
eral sources of uncertainty affect the determination of 
defui- the error on the observed HI mass, the statistical 
uncertainty on the calibrated coefficients and the deter- 
mination of the optical diameter. It is hard to quantify 
a value for the final uncertainty, but a comparison of 
the values derived for a single galaxy using different data 
available in the literature suggests a mean dispersion of 
0.2 — 0.3 dex on the Hl-deficiency. 

Figure [T]a, shows the distribution of the Hl-deficiencies 
for our sample. If we arbitrarily consider unperturbed 
those galaxies which have lost less than a factor of 2 
in their HI mass (de/ni < 0.3), our subsample of iso- 
lated galaxies is composed of 25 objects with a mean 
deficiency defm = —0.045 ±0.27. This is consistent with 
zero. Moreover, we point out that the dispersion on the 
mean value is within the deficiency threshold of 0.3 dex. 



The value for the reference mass Mfja.cxp is computed 
from the observe d correlation betwe en the H2 and the H- 
band luminosity (jBoselli et al.lll99"7l ) for isolated galaxies 
{defm < 0.3). This definition is based on the natural 
scaling of the molecular gas abundance with the size of 
a galaxy, i.e. more massive spiral galaxies are more gas 
rich. Although our sample covers a small range in lumi- 
nosity, owing to the good quality of our data (see section 
14. ip . we choose to compute a new reference mass with a 
least square fit: 



\0gMu2 



0.501ogLH-l-3.98, 



(10) 



Since the dispersion (rms) in the calibration is ^ 0.2 
(See also Fig. [1)3), galaxies with defn^ < 0.2 cannot 
be distinguished from normal galaxies. In the following 
analysis we therefore consider a galaxy H2-deficient only 
if defile ^ 0.2. However, as we shall see below, de/nj 
shows a continuous rather than a bimodal distribution. 
Thus, this criterion is arbitrary. It does have the virtue, 
though, of nicely dividing our sample into classes which 
have a comparable number of objects. As with defm, it is 
hard to quantify an uncertainty for de/na • Deficiency is 
a relative quantity, so the level of uncertainty for defn^ 
is likely to be significantly lower than for the absolute 
H2 mass. The error of 0.7 dex we quoted in equation 
([2|) is therefore an upper limit on the scatter in defn^ . 
A more reliable estimate of the error is the dispersion 
of 0.2 dex we measure in de/na for the control sample 
of isolated galaxies (Fig. [Tja). Since random error in 
computing Mn^, for example galaxy-to-galaxy scatter 
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TABLE 2 

COMPLETNESS FOR THE SURFACE DENSITY PROFILES. 



Fig. 2. — Comparison between calculated defn^ and defm values 
for our sample. The vertical dashed line corresponds to defm = 0.3 
dex, the value below which a galaxy is considered unperturbed; 
the horizontal line represents defu^ = 0.2 dex. The absence of 
galaxies in the upper-left quadrant suggests that Hl-deficiency is 
only a necessary condition for H2-deficiency. 
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Fig. 3. — Average H2 profiles in bins of normalized radius. The 
dashed line is computed by excluding HI— rich galaxies with defm < 
—0.3. The offset among the profiles reflects the defu^ distribution 
in Fig. [2] 



in X, can only broaden the intrinsic dispersion of H2 
masses at fixed L// , our measured scatter of 0.2 dex sets 
a more realistic upper limit on the size of random errors. 
The one caveat is that this value may neglect systematic 
uncertainties if unperturbed and perturbed galaxies have 
different values of luminosity or metallicity, for example, 
and thus systematically different values of X. However, 
we do not observe or expect any obvious trend of these 
quantities as a function of defm or def^^ . 



Class 


Objects 


% Class 


% Total b 


Symbol 


Non-deficient 


10 


40% 


21% 


Circle 


Non gas-rich 


8 


38% 


17% 


Line 


Hl-deficient 


5 


38% 


11% 


Square 


Gas-deficient 


5 


56% 


11% 


Triangle 



^ % Class is defined as the number of galaxies with HI maps 
in each class over the total number of objects in that class. 
^ % Total is defined as the number of galaxies with HI maps in 
each class over 47, the number of objects in the entire sample. 

Symbols used for each class in the plots throughout the pa- 
per. 

Figure [2] shows a comparison between de/na and defm 
in our sample. The vertical dashed line corresponds to 
defm = 0.3 dex, the value below which a galaxy is 
considered unperturbed; the horizontal line represents 
rfe/H2 = 0.2 dex. Inspecting this plot, the first strik- 
ing result is that 9/22 of the galaxies with defm > 0.3 
have defu^ values greater than 0.2 dex, while the re- 
gion characterised by defm < 0.3 and def^^ > 0.2 is 
almost empty. This distribution proves the existence of 
the molecular deficiency in a significant number of galax- 
ies which are HI poor. This figure indicates that defn^ 
and defm do not follow a tight correlation, but this is not 
surprising. In fact, we do not expect effects that deplete 
the HI to act in the same way on the H2 gas which is 
more strongly bound. Rather than a correlation, it is ev- 
ident that galaxies lie in three regions roughly defined 
by our criterion of isolation and molecular deficiency 
(dashed lines). Based on Fig. [21 we divide our sam- 
ple in three classes: unperturbed or non-deficient galax- 
ies {defm < 0.3 and de/na < 0.2), gas-deficient galax- 
ies (defm > 0.3 and de/na > 0.2) and galaxies which 
are only Hl-deficient {defm > 0.3 and def^^ < 0.2). 
This classification should not be considered strict, since 
galaxies are distributed continuously rather than in dis- 
tinct populations, but it is useful for the analysis. We 
test the robustness of our result by changing the isolation 
criterion in the range defm = 0.1 — 0.4 and recomputing 
de/jja after that a new calibration has been derived. We 
find that the distribution in Figure [2] does not change 
significantly. 

Since the definition of the molecular gas deficiency re- 
lies on Eq. (fTO|) which sets the reference mass, the be- 
haviour shown in Fig. [2] could be an artifact imposed 
by this calibration. We further explore the robustness of 
a H2-deficient population by studying the H2-deficiency 
with surface density gas profiles, which is independent 
from the calibration of logMna.exp- In Figure [3] we show 
the average molecular gas surface density profiles, in bins 
of normalized galactocentric radius^ (i?/i?25)- Not all the 
galaxies in our sample are imaged in HI. Therefore, for 
this kind of analysis, we consider only 20 galaxies (See 
Table [1]). Details on the completeness in each class are 
provided in Table [2l By comparing the values for the 

^ This quantity is well defined under the hypothesis that the 
stellar disk is not truncated during the perturbation. Clearly this 
is not the case for tidal interactions. In our sample, NGC 3351 is 
the galaxy most likely to be affected by tidal stripping. However, 
we note that a value for R25 larger than what we have adopted 
would shift the surface density profile at lower H/-R25, enhancing 
the separation between unperturbed and gas-deficient galaxies. 
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Fig. 4. — Far infrared luminosity at 100/.im as a function of the 
molecular gas deficiency. Gas— deficient galaxies show the lowest 
100fj,m luminosity. 

three classes of galaxies (non-deficient with circles, gas- 
deficient with triangles and Hl-deficient with squares), 
it is evident that gas-deficient galaxies have the lowest 
molecular gas column densities, while Hl-deficient and 
non-deficient galaxies exhibit a similar behaviour. The 
dashed line in Fig. [3] shows the average profile excluding 
those galaxies with de/ni < —0.3. The similarity of this 
profile to that of the full non-deficient sample (circles) 
indicates that the observed difference in H2 content be- 
tween the three galaxy classes is not an artifact enhanced 
by the presence of gas-rich galaxies in our sample. Based 
on an independent method. Fig. [3] confirms the existence 
of H2-deficient galaxies in ~ 40% of the HI poor galaxies 
and ~ 20% of all the galaxies in our sample. As discussed 
in the next section, these values should be regarded as 
lower limits. 

4. DISCUSSION 

4.1. Comparison with previous studies 

The first issue that we would like to address is why 
previous studies, based on large CO surveys, did not 
reveal a population of H2-deficient galaxi es. As dis- 
cussed in the i ntroducti on, several s tudies (iStark et al.l 
ToM IKenney fc Youn^TIall flOSl ICasoh et alJ 119911 
Boselli et al.l Il997i . i2Q02 ) based on large samples both 
in the Virgo and in the Coma clusters did not show any 
evidence for a molecular gas deficiency. As discussed by 
[Boselli fc Gavazzi (,20061 . some biases influence the pre- 
vious works; for example, a constant X conversion factor 
was often assumed or samples were selected according to 
FIR selection criteria which prefer CO rich galaxies. Of 
course, our sample has also been selected to have IRAS 
fluxes from the Point Source Catalogue higher than 10 
Jy at 100 /xm. By comparing the lOO^m flux with our 
defu^ in Figured we notice that gas-deficient galaxies 
show systematically lower lOO^m luminosity; therefore, 
we might miss a significant fraction of molecular poor 
galaxies. 



Based on this work, we can identify two additional rea- 
sons for why previous studies failed to detect the molec- 
ular deficiency: the completeness in the flux sampling 
and the homogeneity of the data. While observing at a 
single-dish, the integration time required to detect the 
CO flux is very high (120 minutes at a 12m telescope 
are needed to reach a sensitivity of 3 mK at a resolu- 
tion of 15 km/s) and several observations are required 
to image the entire galaxies. It has only recently be- 
come possible produce CO images in a reasonable ob- 
serving time, thanks to the the availability of multibeam 
receivers that reduce the number of pointings required 
to fully map the entire disk. Due to these technical lim- 
itations, previous studies were based on the observation 
of one position per galaxy plus one or two positions off- 
center, along the major axis. The total CO flux was then 
recovered by interpolating with an exponential profile, a 
procedure that inevitably gives only a rough estimate of 
the prop er H2 mass. By c o mpar ing the integrated CO 
flux from lKennev fc Youngj (|1988D and GOLDMine with 
the Nobeyama and BIMA data, we flnd a dispersion of 
up to 1 dex. In contrast, by comparing galaxies in com- 
mon in the BIMA and Nobayama surveys, the observed 
discrepancy in the total mass is less than 0.4 dex. There- 
fore, we believe that the data used in previous studies 
suffer from fluctuations larger than the observable H2- 
deficiency (< 0.8 dex in our sample). Another possible 
problem is related to the homogeneity of the sample. To 
increase the statistics, earlier studies often combined CO 
data from surveys with different sensitivities and diverse 
techniques to recover the integrated flux, increasing the 
relative dispersion in the observed H2 mass. We attribute 
the fact that the H2-deficiency has remained hidden in 
previous analysis to the combination of these two effects. 
On the contrary, based on high quality and homogeneous 
data from the BIMA and Nobeyama surveys, our work 
allows for a determination of molecular gas deficiency. 

4.2. The observed distribution in the H2-deficiency 

The observed distribution in defn^ deserves a detailed 
analysis. In particular, we are interested in understand- 
ing why Hl-deficiency is a necessary but insufficient con- 
dition for H2-deficiency and what the fundamental dif- 
ference is between galaxies which are only Hl-deficient 
and galaxies which are also molecular gas deficient. The 
answers to these questions lie in the physics that governs 
the H2 formation. Molecular hydrogen is synthesized in 
spiral disks through the interaction of atomic hydrogen 
with dust grains, whereas it is destroyed by the UV radi- 
ation field. The balance between these processes deter- 
mines the ratio of molecula r to atomic hydrogen. In a 
series of recent papers, Kru mholz et al.l pOOB, 2009) in- 
vestigate the problem of atomic to molecular ratio with 
a model that includes the treatment of H2 formation, 
its self-shielding and the dust shielding. According to 
their prescription, the molecular fraction in a galaxy is 
determined primarily by the product of its hydrogen col- 
umn density and metallicity. Under the assumption of a 
two-phase medium (Wolfire et al.l [19951 ) . the results are 
largely independent on the UV flux because a change 
in the UV flux produces a compensating change in the 
HI density. Hence, the effects of UV radiation field are 
largely cancelled out. For the purposes of our analysis, 
the most important prediction of the model by Krumholz 
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Fig. 5. — Average total gas surface density profiles in bins of 
normalized radius, divided according to their de/jji and defu^ 
values. The dashed line is computed excluding HI rich galaxies 
with defiii < —0.3 and the open squares indicate the profile for 
Hl-deficient galaxies without NGC 4321. Gas— deficient galaxies 
have in the disk a gas column density below the values required to 
form molecular hydrogen (~ 10 Mq pc~^). 
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Fig. 6. — Average HI surface density profiles in bins of normal- 
ized radius, divided according to their de/ni and defH2 values. 
The dashed line is computed excluding the Hl-rich galaxies with 
d-efm < —0.3 and the open squares indicates the profile for HI— 
deficient galaxies without NGC 4321. Hl-deficient galaxies have 
in the inner part of the disk an HI column density large enough to 
form H2. 



et al. is that the fraction of the gas in molecular form in- 
creases as the total gas column density S^as = 5]hi +5]h2 
and the metallicity Z increase, providing a large enough 
shield against the photodissociation. In agreement with 
what is observed in local galaxies, the model predicts that 
the HI column density is expected to saturate at a value 
of ^ 10 M0 pc~^. Based on the works by Krumholz et 
al., it is natural to study the behaviour of the total gas 
column density in our sample to understand the nature 
of the Il2-deficiency. 

Figure [5] shows the average total gas surface density 
profiles Shi+Sh2 bins of the normalized galactoccntric 
radius. Inside the optical radius, the total gas content of 
the Hl-deficient galaxies (squares) is the same as that 
of the unperturbed galaxies (circles), while gas-deficient 
galaxies (triangles) have a lower total gas column density. 
This behaviour reflects the observed atomic and molecu- 
lar gas content seen in Fig. [21 yet add new insights into 
the nature of the H2 deficiency. Gas-deflcient galaxies 
have a gas column density too low to shield molecules 
from photodissociation, while Hl-deficient galaxies have 
normal gas content inside their optical radii, which al- 
lows the production of molecular hydrogen as in unper- 
turbed galaxies. Indeed, gas-deficient galaxies lie below 
the critical shielding column density of ~ 10 Mq pc^ 
while unperturbed and Hl-deficient galaxies lie above 
this threshold. As in Fig. [31 we also plot the average 
gas profile excluding galaxies with defm < —0.3 (dashed 
line) : the consistency of the profiles with or without HI- 
rich galaxies excludes any bias in the mean value towards 
high column density. Inspecting the individual HI pro- 
files, we noticed that NGC 4321 exhibits an unusual high 
column density. Excluding this object from the analysis 
(open squares) does not change the results beacause the 
major contribution to the total gas column density comes 



from the molecular hydrogen. 

At this point we have offered an interpretation for why 
the Hl-deficiency is only a necessary condition for the 
in H2-deficiency in HI poor galaxies, but a new question 
arises: if unperturbed and Hl-deficient galaxies have a 
comparable gas column density, why do they have dif- 
ferent values of defm? The analysis presented in Fig- 
ure [S] provides the answer. In this plot we show the 
average HI surface density profiles in bins of the nor- 
malized galactoccntric radius up to twice i?25- As for 
Fig. [51 inspecting the individual HI profiles, we noticed 
that NGC 4321 exhibits an unusual high column den- 
sity. To test if this is relevant for our analysis, we over- 
plot with open squares a new HI profile excluding this 
galaxy. This profile is no longer fully compatible with 
the one observed for non-deficient galaxies. It is also 
evident that the surface density profiles for Hl-deficient 
and gas-deficient galaxies exhibit a different behaviour. 
The former rises towards the center, approaching the gas 
column density observed in unperturbed galaxies, while 
the latter remains at the lowest observed column density 
at all the radii. Assuming the profile for unperturbed 
galaxies serves as a template for how the atomic gas is 
distributed in spiral galaxies, Hl-deficient galaxies have 
a significant residual quantity of gas towards the center 
that allows the H2 formation, but have much less HI at 
larger radii. Therefore, the main contribution to de/jji 
comes from the outer part of the disk where the HI pro- 
file rapidly declines. Gas-deficient galaxies have an HI 
column density that is reduced relative to that of unper- 
turbed galaxies in both their inner and outer disks and 
it is only this reduction in the HI content from the inner 
disk that gives rise to H2 deficiency. The dashed line in 
Fig. [6] shows that the observed behaviour is not driven 
by HI rich galaxies. 
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Fig. 7. — Ratio of the molecular to atomic hydrogen as a function 
of the total gas column density. The solid lines show the expected 
ratio according to the model by Krumholz et al. (2009) for lower 
{[O/H] = —0.4) and upper {[O/H] = 0) value of metallicity in our 
sample. At high gas surface densities gas rich and gas poor galaxies 
exhibit the same behaviour because locally the physics of the H2 
formation remains unchanged. 

In Figure [7] we show the ratio of the molecular to 
atomic hydrogen as a function of the total gas column 
density. In this plot, all systematic differences between 
the classes of galaxy disappear. Since we are explor- 
ing the high density regions, the fact that even gas- 
deficient galaxies show the same behaviour observed in 
non-perturbed and Hl-deficient galaxies suggests that 
molecule formation is a local process within a galaxy and 
that the relevant local variable is column density. The 
solid lines in Fig. [7] are the expected values of mol ecular 
fraction from the model bv lKrumholz et ahl (|2009f ). The 
model is dependent on the metallicity, so we consider a 
lower {[O/H] = -0.4) and an upper {[O/H] = 0) value 
that reasonably correspond to the range of metallicity 
observed in our sample. There are some approximations 
that limit the robustness of the comparison between the 
data and the model. First, the density considered in 
the model is computed assuming pressure balance, which 
may not hold in perturbed galaxies; second, we average 
the gas column density on physical scales that contain 
several atomic-molecular complexes, while the model is 
formulated for a single cloud. Despite all these limita- 
tions, there is a satisfactory agreement between data and 
the model. Finally, Figure[7|is useful to stress once again 
the main result of our analysis: galaxies which suffer from 
HI depletion inside the optical disk are on average shifted 
at lower column density where the production of H2 be- 
comes inefficient. 

4.3. H2-deficiency and star formation 

Our analysis offers a possible solution to the cur- 
rent inconsistency in studies of the star formation and 
galaxy environment. In fact, the observed molecular de- 
ficiency is the missing link between the quenching of 
the star formation rate (SFR) in Hl-deficient galaxies 



(e.g. iGavazzi et al.l l2002f ) and the fact that the bulk 
of the stars form in molecular regions (iWong fc Blitz! 
2002t iK ennicut t et al.ll2007t iFumagalh fc Gavazzill2008r 
Bigiel e t al. 200^. Our simple picture is that, when per- 
turbations reduce the atomic hydrogen column density 
below the threshold required to produce the molecular 
hydrogen, star formation will be suppressed. 

We can test this hypothesis by comparing in Fig- 
ure [8] the integrated atomic and molecular hydrogen 
content with our indicator for the star formation ac- 
tivity. Th e top-left panel rep roduces the well-known 
trend (e.g. IGavazzi et al.l 120021 ) between Hl-deficiency 
and the quenching of star formation activity, measured 
through the Ha + [Nil] E.W.. Gas-deficient galaxies 
(triangles) clearly exhibit a lower specific star formation 
rate than unperturbed (circles) and Hl-deficient galaxies 
(squares). In the top-right panel, we plot the Ha+ [Nil] 
E.W. versus defii^. Two distinct behaviours can now be 
identified: molecular hydrogen deficient galaxies (trian- 
gles) have their star formation rate strongly suppressed 
with respect to molecular rich galaxies (squares and cir- 
cles). Comparing both the top panels, we conclude that 
the trend between de/jji and the star formation rate is in 
reality a reflection of the link between defn^ and de/jji 
with defile and the star formation rate. In more detail, 
looking at the galaxies with defm > 0.3 it appears that 
the Hl-deficient galaxies can have star formation rates 
equal to the values observed in non-deficient galaxies, 
while gas-deficient galaxies have the lowest Ha + [Nil] 
E.W. The dispersion observed in the top-left panel can 
be interpreted as being primarily a product of the ob- 
served range in the molecular gas deficiency: not all Hl- 
deficient galaxies are molecular deficient and thus not all 
them have suppressed star formation. This degeneracy 
is removed by directly comparing the Ha + [Nil] E.W. 
with the rfe/fja • Because we find the same result by com- 
paring the 100/zm luminosity normalized to the H-band 
luminosity with both defm and defn^ (bottom panels), 
we are confident that dust extinction or [Nil] contami- 
nation do not significantly affect our conclusions. 

5. SUMMARY AND CONCLUSION 

Using high sensitivity and high resolution CO images 
from the Nobeyama and the BIMA surveys we have stud- 
ied molecular hydrogen deficiency in spiral galaxies. Our 
results and conclusions are summarised as follows. 

- Comparing the total masses of molecular and 
atomic hydrogen, we detect H2-deficiency in a sub- 
set of galaxies which are HI poor. The same result 
is found by studying directly the H2 surface density 
profiles. We observe that molecular gas deficiency 
is associated with the HI reduction inside the op- 
tical disk. The deficiency in the atomic hydrogen 
is a necessary condition for the molecular gas defi- 
ciency, but the de/ua parameter spans a wide range 
of values even in galaxies which arc Hl-deficicnt. 
In fact, there are HI poor galaxies which have a 
normal H2 content. 

- When we examine the molecular fraction at a given 
point in a galaxy as a function of the total gas 
column density, we no longer detect any system- 
atic difference between unperturbed, Hl-deficient 
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Fig. 8. — Comparison between the atomic and molecular gas deficiency with SFR indicators. The top panels show the Ha + [Nil] E.W., 
while the bottom panels show the 100/^m luminosity normalised to the H-band luminosity. Gas-deficient galaxies have the lowest absolute 
specific star formation rate as expected if stars form in molecular rich regions. 

and gas-deficient galaxies. This suggests that the 
total gas column density, or something strongly 
correlated with it, is the primary factor in deter- 
mining the molecular content of a galaxy. In this 
picture, a galaxy moves from Hl-deficicnt to gas- 
deficient when it loses enough gas inside its opti- 
cal disk to bring the gas column density below the 
critical value of ~ 10 Mq pc~^ required for H2 to 
be present. This behaviou r is quantitative l y con- 
sisten t with the models of iKrumholz et aD ()2008l . 
12009^ . albeit with significant uncertainties. 

- Since the bulk of the star formation seems con- 
nected with the molecular gas, when perturbations 
reduce the atomic hydrogen column density below 
the threshold required to produce molecular hydro- 
gen, star formation is suppressed. Studying the 
Ha + [Nil] E.W. or the 100/zm luminosity, we 
show that H2-deficient galaxies have less star for- 
mation than what is observed in H2-normal galax- 
ies and that the well-known trend between the HI- 
deficiency and the star formation rate most likely 
reflects a physical connection between molecular 
hydrogen depletion and star formation activity. 

All together, these results corroborate a picture in 
which the gas phase as a whole feeds the star forma- 
tion activity. The atomic hydrogen is the essential and 
primary fuel for the formation of molecular hydrogen 
that sustains the bulk of star formation activity. When 
environmental processes or feedback are able to reduce 
the atomic gas inside the optical disk, the HI column 
density drops below the critical column density required 
for molecular hydrogen formation and star formation is 



quenched. This scenario is supported on cosmological 
scales, where the HI content in gal axies appears to be 
unchanged for the past ~ 10 Gyr (jProchaska fc Wol^ 
I2OO8I ). This suggests that the global star formation rate 
is driven by the accretion of fresh gas from the IGM. In 
addition, our analysis explains how environmental per- 
turbations or galaxy feedback, which remove the atomic 
gas, can suppress the formation of the H2 and stars. The 
behaviour seen in the innermost bin (0 < R/R25 < 0.2) 
of Figures [H [5] and [6] allows to rule out the hypothesis 
that AGN feedback is the dominant process responsi- 
ble for the gas depletion. In fact, gas is almost unper- 
turbed in the galaxy center. This is expected in our 
sample where 7/9 of gas-deficient galaxies belong to the 
Virgo cluster or the Coma 1 cl oud in which gas strip- 
ping is the most eff icient process (|Boselli fc Gavazzill2006l : 
IBoselh et alll2006h . 

Ongoing CO surveys as HERACLES by Leroy et al. 
(in prep) or JCMT, th e Nearby Galaxies Legacy Survey 
bv lWilson etld] ((2008h . will soon provide larger samples 
to test our results. A preliminar y comp arison with the 
fluxes published by iWilson et al.l (|2008f ) appear to sup- 
port our work. Looking at their Table 4, the four listed 
galaxies show from the top to the bottom decreasing val- 
ues of SFR and CO J : 3 — 2 luminosity; consistently, 
we flnd for them an increasing value of defn^. If con- 
flrmed, our results support the argument that the bulk 
of star formation is connected with the molecular gas 
phase. Moreover, the ability to derive the molecular gas 
content starting from an observed HI distribution will 
provide a useful way to correct the sta r formation rate 
expected from loc al star formation laws (|Kennicuttlll998l : 
iBigiel et al.ll2008f l. including the effect of the molecular 
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gas deficiency which suppresses star formation in rich 
clusters or in interacting galaxies. 

Our analysis leaves open a fundamental question on the 
implications that H2-deficiency might have for galaxy 
evolution. In the standard scenario where star for- 
mation is supported by cold gas infall towards the in- 
ner part of the galaxy halo (jBoissier fc Prantzoi 119991 : 
iDekel fc Birnboimll2006[ ). it is fascinating to think that, 
when the gas flow is halted, the gas is devoured by on- 
going star formation. The galaxy depletes its gaseous 
reser voirs and turns into a passive and quiescent sys- 
tem (|van den Berghl 119761 ). This process might explain 
a first step for the transition between the blue and red 
sequence in the color-magnitude diagram. Interactions 
with the intergalactic medium might thus be an alter- 
native process to AGN feedback or galaxy mergers to 
explain the departure of star forming systems from the 
blue sequence. In our scenario, two time-scales govern 
this migration, one set by the atomic gas depletion and 
the other related to the response of the star formation ac- 
tivity to the gas removal. The first varies according to the 
physical process responsible for the HI gas removal, but 
in this phase the galaxy is still forming stars (e.g. Boselli 
et al 2006). The second is a relatively short time-scale, at 
most set by the gas co nsumption through star formation 
(- 2 Gvr. lBigiel eraI1 [2008i^ but more likely connected 
with the much shorter life time of the GMCs. This sec- 
ond phase is consistent with a rapid evolutionary path in 
the color-magnitude diagram. However, although it has 
been proved that dwarf galaxies can turn int o dwarf ellip- 
ticals within ~ 2 Gyr in the local universe (|Boselli et al.l 
[200l . it is challenging to prove that massive spirals can 

® This value comes from the ratio of gas surface density and star 
formation rate surface density which is assumed constant. How- 
ever, gas recychng may increase Egas up to 40% and the actual 



turn into giant ellipticals just as a consequence of gas 
removal through interactions with the cluster IGM. In 
fact, their angular momentum would be fairly conserved 
in this kind of interaction. Major mergers a s those dom- 
inating at higher redshifts must be invoked (jFaber et al.l 
,2007i) . 
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time-scale may be longer. 



APPENDIX 
HI FLUXES 

As discussed in Sec. 12.2.21 the HI fluxes quoted in the literature for individual galaxies are subject to discrepancies 
which in some case are not negligible. To show that our result is independent of a particular cho ice for the HI fluxes, w e 
compute different sets of de/ni parameters using fluxes from ALFALFA, THINGS, VIVA and lSpringob et all (pOOSh . 
Figure [9] shows a comparison of these values with the Hl-deficiency computed from HyperLeda fiuxes adopted in our 
analysis. 

Inspecting the left-panels, we notice that defm computed from ALFALFA and VIVA is in agreement with our 
values. Because these two surveys provide HI fluxes in the Virgo cluster where the most Hl-deficient objects reside, 
we conclude that defm for these galaxies is not affected by important fluctuation. On the contrary, looking at the 
right -panels, we notice that at low Hl-deficiencies a systematic offset appears. Since only few galaxies at low deficiency 
are available in ALFALFA or VIVA, an additional comparison is almost impossible and it is difficult to further discuss 
this difference. However, assuming that HyperLeda fluxes suffer from a systematic offset, we should correct our defm 
towards even lower defm values and galaxies we classify as non-deficient remain in this class. This is not the case 
for NGC 5236 and NGC 6946, the two outliers visible in the bottom-right panel. However, these galaxies are more 
extended than the primary beam mapped by THINGS and we are confident that our value for defm is more accurate. 

This comparison conflrms that, despite the fluctuation in the defm parameter, the flnal distribution for the HI- 
deflciency and thus our analysis are not significantly affected by a particular choice of the HI flux. 
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